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ABSTRACT 

Context. We consider a localized pulse in the component of velocity, parallel to the ambient magnetic field lines, that 
is initially launched in the solar chromosphere. 

Aims. We aim to generalize the recent numerical model of spicule formation (Murawski & Zaqarashvili 2010) by 
implementing a VAL-C model of solar temperature. 

Methods. With the use of the code FLASH we solve two-dimensional ideal magnetohydrodynamic equations numerically 
to simulate the solar macrospicules. 

Results. Our numerical results reveal that the pulse located below the transition region triggers plasma perturbations, 
which exhibit many features of macrospicules. We also present an observational (SDO/AIA 304 A) case study of the 
macrospicule that approximately mimics the numerical simulations. 

Conclusions. In the frame of the model we devised, the solar macrospicules can be triggered by velocity pulses launched 
from the chromosphere. 
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1. Introduction 

Spicules are thin, cool and dense structures that are ob- 
served in the solar limb (Beckers I1968| I1972[ Suematsu 
H9981 Sterling [20001 Zaqarashvili & Erdelyi [2009) . Type I 
spicules are seen to be formed at an altitude of about 2000 
km where they reveal a speed of 25 km s -1 , reach at a 
maximum level and then either disappear or fall off to the 
photosphere. A typical lifetime of spicules is 5 — 15 min with 
an average value of ~ 7 min (PasachofT et al. 2009). Spicules 
consist of double thread structures (Tanaka 1974, Dara et 
al. I1998[ Suematsu et al. l2008j) and reveal the bi-directional 
flow (Tsiropoul a et a l. H9941 Tziotziou et al. [20031 [20041 
PasachofT et al. 2009). Typical electron temperature and 
electron density in spicules are respectively (15 — 17) • 10 3 
K and 2-10 11 - 3.5 • 10 10 cm" 3 at a ltitudes of 4-10 Mm 
above the solar surface (Beckers |1968|) . As a result, spicules 
are much cooler and denser than ambient coronal plasma. 
High resolution observations by Solar Optical Telescope on- 
board Hinode revealed another type of spicules with many 
features different than classical limb spicules and they are 
referred as type II spicules (De Pontieu et al. 2007). Mean 
diameter of type I spicules is estimated as 660 ± 200 km 
(PasachofT et al. l2009p . but the type II spicules have smaller 
diameters (< 200 km) in the Ca II H line (De Pontieu et al. 
2007). While observed in the R a line, type I spicules can 
reach up to 4 — 12 Mm in height from the solar limb with 
a mean value of 7200 ± 2000 km (PasachofT et al. 12009]) . 
The type II spicules are significantly shorter. The main dif- 
ferences between type I and type II spicules are: (a) the 



lifetime, the former lasts longer; (b) the evolution, the for- 
mer usually shows a parabolic profile in Ca II H while the 
latter shows an upflow and then disappears; (c) the for- 
mer usually reveal smaller velocities (e.g., De Pontieu et al. 
[20071120091 Rouppe Van der Voort et al. l2009]) . 

Additionally, very long spicules, called as macrospicules 
with typical length of up to 40 Mm and with higher tem- 
perature are frequently observed mostly near the polar re- 
gions (Cook et al. 119841 Pike & Harrison 119971 Ashbourn 
& Woods 12005]) . Georgakilas et al. ([T999]) have observed 
the giant solar spicules that reach at a maximum height 
of about 12 — 15 Mm with their lifetime of the order of 
5 — 12 minutes, and ejection velocity around ~50 km s _1 . 
Wilhelm (2000) has also observed the largest macrospicules 
upto a height of 40 Mm in the north polar coronal hole. In 
conclusions, the solar macrospicules can extend from 7 Mm 
to 45 Mm above the solar limb with lifetime of 3-45 min, 
and they are mostly concentrated in the polar coronal holes 
(Sterling [2000]). 

In spite of various theoretical models to explain the 
spicule ejection in the lower solar atmosphere, many recent 
numerical methods have been developed to simulate the 
solar spicules/macrospicules with an energy input at their 
base in the photosphere as a pressure pulse or an Alfven 
wave that steepens into a shock wave (Sterling 20®] and 
references therein). Recently, Hansteen et al. (2006J and 
De Pontieu et al. (2007) simulated the formation of dy- 
namic fibrils due to slow magneto-acoustic shocks through 
two-dimensional (2D) numerical simulations. They suggest 
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that these shocks are formed when acoustic waves gener- 
ated by convective flows and global p-modes in the lower 
lying photosphere leak upward into the magnetized chromo- 
sphere. Heggland et al. (j2QQ7|) used the initial periodic pis- 
ton to drive the upward propagating shocks in ID simula- 
tion and Martinez- Sykor a et al. (2009) considered the emer- 
gence of new magnetic flux but the drivers of spicules come 
from collapsing granules, energy release in the photosphere 
or lower chromosphere. However, these simulations could 
not mimic the double structures and bi-directional flows 
in spicules. On the other hand, Murawski & Zaqarashvili 
(2010) have performed 2D numerical simulations of mag- 
netohydrodynamic (MHD) equations and showed that the 
2D rebound shock model (Hollweg 1982) may explain both 
the double structures and bi-directional flows. They used 
a single initial velocity pulse, which leads to the formation 
of consecutive shocks due to nonlinear wake in the strati- 
fied atmosphere. However, they considered a simple model 
of atmospheric temperature which was approximated by a 
smoothed step function profile. 

In this paper, we adopt a more appropriate tempera- 
ture profile of Vernazza et al. (1981) that implements char- 
acteristic (Alfven and tube) speed profiles and launch the 
pulse at the chromosphere. Our more general model re- 
veals that macrospicules are effectively excited by veloc- 
ity pulses launched at the chromosphere. Additionally, our 
model confirms the findings of Murawski & Zaqarashvili 
(|2010|) who explained the observed properties of spicules 
(e.g., width, height, and bi-directional flows). We also 
present an observational case study using SDO/AIA 304 
A of a macrospicule in the north polar coronal hole 
(NPCH), which approximately matches with the simu- 
lated macrospicule in oblique magnetic field. This observed 
spicule appears near the north-east boundary of a moder- 
ately evolved hole at the polar cap. 

This paper is organized as follows. A numerical model 
is presented in Sect. 2, and the corresponding numerical 
results are shown in Sect. 3. The observational case study 
of a macrospicule is presented in Sect. 4 in support of the 
numerical simulations. This paper is completed with a dis- 
cussion and conclusions in Sect. 5. 



2.1. Equilibrium configuration 

We assume that the solar atmosphere is in static equilib- 
rium (V e = 0) with a force-free magnetic field, 

(VxB e ) xB e = 0. (5) 

As a result, the pressure gradient is balanced by the gravity 
force, 

- Vp e + Qeg = . (6) 

Here the subscript e corresponds to equilibrium quantities. 
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2. A numerical model 

We consider a gravitationally-stratified solar atmosphere 
which is described by the ideal 2D MHD equations: 



| + v.(,v) = o, 
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-Vp+ -(V x B) x B + £g, 

g + V.(pV) = (l-7)pV-V, P=^QT, 
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— = Vx(VxB), V-B = 0. 
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Here q is mass density, V = [V^,V^,0] is flow velocity, 
B = [£> x , By, 0] is the magnetic field, p is gas pressure, T is 
temperature, 7 = 5/3 is the adiabatic index, g = (0, —g, 0) 
is gravitational acceleration of its value g = 274 m s -2 , m 
is mean particle mass and &b is the Boltzmann's constant. 



Fig. 1. Equilibrium profile of solar temperature (top panel) and 
mass density (bottom panel). 



Using the ideal gas law and the ^-component of hydro- 
static pressure balance indicated by Eq. |6|, we express 
equilibrium gas pressure and mass density as 



Pe(y) =Po exp 
Here 

A(y) = k B T e (y)/(mg) 



dy 

W) 



Qe(y) 



Pe(y) 
9Mv) 



(7) 



(8) 



is the pressure scale- height, and po denotes the gas pressure 
at the reference level that we choose in the solar corona at 
y T = 10 Mm. 

We adopt an equilibrium temperature profile T e (z) for 
the solar atmosphere that is close to the VAL-C atmo- 
spheric model of Vernazza et al. (1981). It is smoothly ex- 
tended into the corona (Fig. [I] top). Then with the use of 
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Fig. 2. The static grid system with block boundaries repre- 
sented by solid lines for the case of xo — 5 Mm. The initial per- 
turbation of the system is displayed by the red dot at (5.00, 1.75) 
Mm. 



For these settings magnetic field is predominantly ver- 
tical around x = Mm, while around [x = 5, y = 1.75) 
Mm it is oblique with about 7r/3 angle to the solar surface. 
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Fig. 6. Same as Fig. 4, but for the case of oblique magnetic 
field shown in Fig. 5. 



Eq. Q we obtain the corresponding gas pressure and mass 
density profiles. 

We assume that the initial magnetic field satisfies a 
current-free condition, V x B e = 0, and it is specified by 
the magnetic flux function, A, such that 



with 

A(x : y) = S 1 



dA dA 

dy ' dx 



,0 



x — Xi 



(x - xi) 2 + (y- yi) 2 ' 



(9) 



(10) 



Here Si denotes a strength of magnetic moment that is 
located at (xi, y±). We choose and hold fixed x\ = Mm, 
yi — —40 Mm, and Si is specified from the requirement 
that at the reference point (x T = 0, y r = 10) Mm Alfven 
speed, c A (x Y ,y r ) = |B e (x r , y v )\/ \J V>Qe(yr), is 10 times larger 
than the sound speed, c s (y r ) = y/ 7p e (2/r) / Qe(Vr) • 
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Fig. 4. Time signature of velocity V y collected at (x = 0, y = 
14) Mm for the case of essentially vertical magnetic field as 
shown in Fig. 3. 



2.2. Perturbations 

We initially perturb the above equilibrium impulsively by a 
Gaussian pulse in the component of velocity that is nearly 
parallel to ambient magnetic field lines, V\\, viz., 



V\\(x,y,t = 0) = A v exp 



(x - x ) 2 + (y - yof 



(11) 



Here A v is the amplitude of the pulse, (#o,2/o) is ^s initial 
position and w denotes its width. We set and hold fixed 
A v = 40 km s _1 , w = 150 km, and yo — 1.75 Mm, while 
allow xq to vary. We consider two cases: (a) xq = Mm 
and (b) xq — 5 Mm. 

3. Results of numerical simulations 

Equations Q-Q are solved numerically using the code 
FLASH (Lee & Deane 2009). This code implem ents a 
second-order unsplit Godunov solver (e.g., Murawski I2002j) 
with various slope limiters and Riemann solvers, as well 
as adaptive mesh refinement (AMR). We use the minmod 
slope limiter and the Roe Riemann solver (e.g., Toro 2009). 
For the case of (a) and (b) we set the simulation box as 
(5, 5) Mm x (1, 41) Mm and (0, 30) Mm x (1, 31) Mm, respec- 
tively. We impose boundary conditions by fixing in time all 
plasma quantities at all four boundaries to there equilib- 
rium values. In all our studies we use a static but non- 
uniform grid with a minimum (maximum) level of refine- 
ment set to 1 (7). See Fig. |5] for the grid system for the case 
of xo = 5 Mm. A similar grid system was chosen for the 
case of xq = Mm. As the grid was static no refinement 
strategy was adopted. As each block consists of 8 x 8 identi- 
cal numerical cells, we reach the effective spatial resolution 
of about 60 km. 

We discuss two cases: (a) essentially vertical magnetic 
field and (b) oblique magnetic field. For the case (a) we 
launch the initial pulse at xq = Mm, while the case (b) 
is realized for xq = 5 Mm. In the latter case magnetic field 
creates the angle of about 7r/3 to the horizontal direction. 
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3.1. Essentially vertical magnetic field 

The temporal snapshots of the simulated macrospicule in 
the model atmosphere of the VAL-C temperature and es- 
sentially vertical magnetic field are shown in Figure [3j It 
displays the spatial profiles of plasma temperature (colour 
maps) and velocity (arrows), resulting from the initial ve- 
locity pulse that was launched at xq = Mm. The upper 
left panel corresponds to t = 100 s. Cold chromospheric 
plasma lags behind the shock front, which at this time 
reaches altitude of y c± 8 Mm. The reason for the material 
being lifted up is the rarefaction of the plasma behind the 
shock front, which leads to low pressure there (Hansteen et 
al. I2006j) . As a result, the pressure gradient works against 
gravity and forces the chromospheric material to penetrate 
the solar corona. 

It is noteworthy that small-amplitude vertically prop- 
agating waves are described by the Klein-Gordon equa- 
tion (Lamb I1909j) . It follows from this equation that ini- 
tial, spatially-localized perturbation results in a wavefront 
which moves away from the launching place. The distur- 
bance ahead of the wavefront is at rest, but behind the 
wavefront an oscillation wake results in. This wake oscil- 
lates with essentially constant frequency and the amplitude 
of the oscillations behind the wavefront declines gradually 
in time. As a consequence of the equilibrium mass den- 
sity fall off with height, the amplitude of these oscillations 
grows with altitude. Then, the linear theory is not valid 
anymore as nonlinear effects result in trailing shocks with 
the secondary shock following the leading shock (Hollweg 
H9821) . 

The next snapshot (top middle panel) is drawn for 
t = 200 s, when the macrospicule raised to the altitude of 
y ~ 11 Mm. At t = 300 s (top-right panel) the macrospicule 
already subsided. As the secondary shock lifts up the chro- 
mospheric material (see two small peaks on both sides of the 
macrospicule), there are upward flows at the macrospicule 
sides, well seen at t = 400 s (middle-left panel). It is clear 
from the snapshots at t > 200 s that the plasma exhibits 
some horizontal oscillations. This is because the plasma 
pushes the magnetic field lines outside due to the large value 
of plasma /3 at the chromosphere. The signatures of these 
horizontal oscillations are the surface waves which are run- 
ning along transition region-corona interface. There are well 
seen bi-directional flows in the macrospicule: downward at 
the centre and upward at the boundaries, which agrees with 
the observational data (Tsiropoula et al. 11994} Tziotziou et 
al. 120031 120041 Pasachoff et al. 12009]) . Eventually, in case of 
our numerical simulation, the downward streaming plasma 
is also following a path in the cool core region as seen at 
t = 800 s. Since the macrospicule plasma at its leading edge 
enters in the corona, it may be heated quickly at the coronal 
temperature, and associated with the greater upflow veloc- 
ity typically in the range of 20 — 40 km s _1 compared to 
the background in case of both singular and bi-directional 
flows (cf., Fig. 3). Therefore, this segment of macrospicule 
may appear in form of Doppler blue-shift of the centroid 
of FUV/EUV lines formed with maximum ionic fraction 
at their typical coronal/TR formation temperatures. The 
blue shift and plasma upflows with the typical macrospicule 
rise-up speed have been observed by Wilhelm (2000]) and 
Suematsu et al. (1995). However, the scenario of temper- 
ature, flow, and magnetic field structures in the bottom 
part of the macrospicules are rather complex. The cool 



core, steep temperature gradient towards boundary may 
evolve at the base. The up and down-flow structuring (cf. 
Fig. |3| may be observed in form of mixed red and blue-shift 
scenario as previously explained by Wilhelm (2000) in the 
pass-band of cool chromospheric lines. 

The middle-right panel is drawn for t = 600 s, which 
clearly resembles a macrospicule-like structure with the 
chromospheric temperature and mass density. Its width and 
mean rising speed are about 0.4 Mm and 30 km s _1 , respec- 
tively. These values are close to the corresponding char- 
acteristics of typical macrospicules. However, this is just 
a presentation of the simulation case study of the solar 
macrospicule in essentially vertical magnetic field config- 
uration. The typical length, width, speed, and life time can 
vary by tuning the initial conditions of the numerical model. 
Afterwards the material again falls back at the central part, 
while the next shock forces the plasma to move upwards at 
the spicule boundaries. We see quasi-recurrently accuring 
single (t = 200 s, t = 600 s, t = 800 s), double (e.g., t = 400 
s), and triple (e.g. t = 300 s) structures which are charac- 
teristic features of macrospicules. 

Figure [4] illustrates the vertical component of velocity 
that is collected in time at the detection point (x = 0, y = 
14) Mm for the case of Fig. [3j As a result of fast mass den- 
sity fall off with height upwardly propagating waves grow 
in their amplitude and steepen rapidly into shocks. The ar- 
rival of the first shock front at the detection point, y = 14 
Mm, is clearly seen at t ~ 70 s. The second shock front 
reaches the detection point at t ~ 460 s, i.e., after ~ 290 
s. This secondary shock results from the nonlinear wake, 
which lags behind the leading shock. 



3.2. Oblique magnetic field 

The temporal snapshots of the simulated macrospicule 
in the model atmosphere of the VAL-C temperature and 
oblique magnetic field are shown in Figure [5] It illustrates 
the spatial profiles of plasma temperature (colour maps) 
and vertical velocity (arrows) resulting from the initial ve- 
locity pulse that was launched at xq = 5 Mm. As a result, 
the plasma flows along the inclined magnetic field. The dy- 
namics of plasma essentially resembles that obtained by 
vertically propagating pulse. At t = 200 s the velocity sig- 
nal reaches the altitude of up to 8 Mm. Then it falls back 
again resembling the consecutive double structure and bi- 
directional flows. 

Figure [6] illustrates the vertical component of veloc- 
ity that is collected in time at the detection point (xd = 
8.5, = 14) Mm for the case of Fig. [5] The arrival of the 
first shock front at the detection point is clearly seen at 
t ~ 100 s. The second shock arrives to the detection point 
at t = 430 s i.e., after ~ 330 s. Another two shocks appear 
at latter times that is at ~ 650 s and ~ 900 s. 

However, this is also a presentation of the simulation 
case study of the solar macrospicule in oblique magnetic 
field configuration driven by a velocity pulse in the model 
temperature of the solar atmosphere. The typical length, 
width, speed, and life time of the obliquely propagating 
macrospicules can also vary by tuning the initial conditions 
of the numerical model. 
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4. Observational evidence of a pulse driven 
macrospicule 

In the present section, we describe an observational signa- 
ture of a limb spicule which approximately resembles the 
simulated macrospicule as driven by a velocity pulse in the 
more general oblique magnetic field configuration. We use 
a time-series data of a solar spicule at north polar coronal 
hole as observed in 304 A filter of Atmospheric Assembly 
Imager (AIA) onboard the Solar Dynamics Observatory 
(SDO) on 3 August 2010 during 03:12:50 UT-03:22:50 UT. 
The SDO/AIA has a typical resolution of 0.6" per pixel 
and the highest cadence of 12 s, and it observes the full 
solar disk in three UV (1600 A, 1700 A, 4500 A) and seven 
EUV (171 A, 193 A, 211 A, 94 A, 304 A, 335 A, 131 A) 
wavelengths. Therefore, it provides the unique observations 
of multi-temperature, high-resolution, and high-temporal 
plasma dynamics all over the Sun. The field-of-view of the 
north polar coronal hole as observed by SDO/AIA 304 A on 
03 August 2010 was (1094", 317"), while the (X cen , Y cen ) 
was (-29.618", 918.631"). The time series has been obtained 
by the SSW cutout service at LMSAL, USA, which is cor- 
rected for the flat-field and spikes. We run aia_prep subrou- 
tine of SSW IDL also for further calibration and cleaning 
of the time series data. 

We observe a giant spicule dynamics at the limb of north 
polar coronal hole. The macrospicule has been launched 
obliquely from the background open field lines of the po- 
lar coronal hole. The life-time of the spicule was observed 
as —10 min, which fits with the typical life time of the 
spicules/macrospicules (Georgakilas et al. 1999, Sterling 
2000, Pasachoff et al. 2009). The observed macrospicule 
reaches upto a height of —12 Mm and has the width of 
—2 Mm. It should be noted that the present observation is 
only a case study in the support of the numerical simulation 
of a solar macrospicule. Macrospicules may reach to 7-40 
Mm heights with a typical life time of 3-45 mins (Sterling 
2000} . Our aim here is only to show that these observations 
reveal upto some extent the nature of the spicule as per the 
numerical modeling. 

The spicule starts moving up above the solar limb on 
—03:13 UT and reaches at a maximum height of —12 Mm 
at —03:17 UT with approximate average rising speed of 
-45 km s" 1 . The height of -12 Mm, life-time of -10 
min, and the speed of —45 km s _1 are typical for the 
observed macrospicules (Georgakilas et al. 1999, Sterling 
2000). The unique evidence this macrospicule presented is 
its rise obliquely off the limb during —250 s of its life time, 
and then its fall back along the same path. The simulated 
macrospicule in the oblique magnetic field configuration 
(cf., Fig. 5) also reaches at a height of —8 Mm in —200 s, 
and then its material is attempting to fall back. The falling 
material of the simulated macrospicule is encountered by 
the wave train of the another upcoming pulses that up- 
lift the spicule material quasi-periodically upto a certain 
height again and again in the solar atmosphere. The ob- 
served dynamics of macrospicule also supports that it may 
be formed by a similar velocity pulse that triggered in the 
chromosphere and steepens in form of a shock in the transi- 
tion region/corona. It should be noted that the rising time 
of the observed macrospicule is —250 s, however, its falling 
time is larger as —350 s. This means that the falling cool 
plasma may be interacting with the uplifted plasma coming 
from the lower solar atmosphere along the same path of the 



spicule. It is also clear from the snapshots during 03:16 UT- 
03:23 UT that the falling plasma is trying to settle down, 
while some brightened material is being pushed up from 
the lower part of atmosphere. However, the signal in the 
line 304 A filter is not necessarily optically thick, and other 
contributions from the line-of-sight signal may contaminate 
the observations. Therefore, it is not completely clear that 
the brightening at the footprint is due to the rebound com- 
pression. After reaching in the upper atmosphere along the 
magnetic field lines, this spicule falls back and cool plasma 
traced back its backward path as observed here. 

Despite of these uncertainties, the observed feature is 
approximately similar to the simulated macrospicule in the 
oblique magnetic field. In the case of the observed spicule, 
finally the plasma is being settled down. The reason may be 
that the upcoming pulses in the real atmosphere may not 
have the same amplitude and spatio-temporal distribution 
at the base of the spicule. The decay of the wave train 
of the pulses does not probably allow the quasi-periodic 
rise and fall of the spicule plasma from the same place. 
The downflowing material may also suppress the upcoming 
pulses that may finally allow the settling down of the spicule 
plasma. 

The observed spicule is two times wider than the sim- 
ulated one. However, the mean rising time and maxi- 
mal height are almost closer for observed and simulated 
macrospicules. The height of the simulated spicule can be 
tuned exactly to the observations by small adjustment of 
the velocity pulse. Simulated and observed velocities are 
also matching to the typical velocities of the macrospicules. 
The fine tuning of the initial conditions in the numerical 
simulation can exactly mimic the observed conditions, how- 
ever, here our aim is only to present an observational case 
study of the macrospicule to reveal its pulse driven nature. 
The resolution and real complex conditions of the solar at- 
mosphere do not permit us exactly to precisely study the 
quasi-periodic rise and fall of the spicule plasma at the same 
place in the observations as were evident in the simulations. 
However, greater life time of the falling plasma, and some 
observational evidence in the time series (03:17-03:23 UT) 
indicate interaction of falling material with the element of 
rising plasma. This may be due to the arrival of another 
pulse from below the chromosphere, which may try to push 
the plasma opposing the gravitational free fall. 

In conclusions, the observations of the macrospicule re- 
sembles the simulation of pulse driven macrospicules. One 
more interesting point is noticeable in the observations that 
the spicule is being generated near the boundary of a bright- 
ened network at north polar coronal hole. This network may 
be the place of the launching of a velocity pulse that trig- 
gers the macrospicule higher in the solar atmosphere. It 
should be noted that all the spicules/macrospicules do not 
show the similar properties as we observe here. Therefore, 
the activity in the associated brightened magnetic network 
may be the most probable cause to drive the initial velocity 
pulse to trigger the observed macrospicule. However, this is 
just a case study to support the numerical simulations pre- 
sented in the paper, and the measured spatial and temporal 
scales are approximate when compared to the simulations. 
The detailed observational search of the solar macrospicules 
matching exactly with the numerical simulations, will be 
done in our future projects. 
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5. Discussion and conclusions 

The formation of solar macrospicules is still an unresolved 
problem in solar physics. Several competitive mechanisms 
have been supposed from time to time, but none of them 
could explain all properties of macrospicules. We performed 
2D numerical simulation of the velocity pulse, which was 
launched at the chromosphere, in stratified solar atmo- 
sphere with the VAL-C temperature profile. The ampli- 
tude of upward propagating perturbation rapidly grows 
with height due to the rapid decrease of equilibrium mass 
density. Therefore, the perturbation quickly steepens into 
the shock in upper regions of the solar chromosphere that 
launches the cool spicule material behind it (Hansteen et 
al. I2006p . The simulated properties of the macrospicules 
both in the vertical and oblique magnetic field configu- 
rations, e.g., velocity, height, life-time, bi-directional flows 
etc, clearly match with the typical properties of such type 
of giant spicules. 

We also present a case study of the observed 
macrospicule using the recent observations from SDO/AIA 
in 304 A to support the observations. The observed 
macrospicule rises obliquely off the limb in north polar 
coronal hole and reaches at a height of ~12 Mm in the 
first 3-4 minutes of its total life time of ~10 min. These ob- 
served phenomena and parameters support the simulated 
macrospicule in the oblique magnetic field configuration. 
The width of the observed spicule is, however, compara- 
tively larger that the simulated one. But, this is one of the 
free parameters of the numerical simulation that can be 
tuned with the shape and size of the velocity pulse to trig- 
ger such macrospicules. The amplitude of the velocity pulse 
can also tune the rising speed of the simulated spicule. The 
observed spicule rises up and then falls back along the same 
path just like the simulated pulse driven macrospicule. The 
doubly splitted head of the observed spicule is also evi- 
dent in the observations once it reaches near the maximum 
height. This may be the indication of the initiation of bi- 
directional plasma flows along the spicule as similar to the 
numerical simulations. Although, the baseline of the obser- 
vations does not allow us to examine the interaction of the 
falling spicule material with the uplifting plasma. However, 
the longer falling phase compared to the rising one may 
indicate some interaction between falling and rising plas- 
mas. Unfortunately, we do not have the Doppler velocity 
fields and related observations to comment more on this 
aspect. However, we can easily identify some upward mo- 
tion of dense bright plasma near the spicule base around 
03:16-03:21 UT that may oppose the falling plasma that is 
already being lifted upto a maximum height of ~12 Mm. 
When such interaction occurs, the falling plasma is proba- 
bly being subsided and shows the two splitted heads as a 
double component as also reported in the simulation snap- 
shot of 400 s in Fig. 5. 

In conclusions, our numerical simulations of the solar 
macrospicules exhibit approximately observed properties. 
We should admit that the performed 2D simulations are 
ideal in the sense that they do not include radiative trans- 
fer, thermal conduction along field lines etc. The field con- 
figuration and the initial stratification are simple and in 
hydrostatic equilibrium. Since we did not attempt to syn- 
thesize observables, the exact comparison with observations 
is unlikely. The detailed observational search of the solar 
macrospicules matching exactly with the numerical simu- 



lations, will be done in our future projects. Simultaneous 
imaging and spectroscopic observations should also be per- 
formed in future to explore extensively the dynamics of such 
macrospicules that will impose rigid constraints on the nu- 
merical simulation models also. 
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Fig. 3. Temporal snapshots of simulated macrospicule in case of essentially vertical magnetic field (see Sect 3.1). The temperature 
(colour maps) and velocity (arrows) profiles at t = 100 s, t = 200 s, t = 300 s, t = 400 s, t = 500 s, t = 600 s, t = 700 s, t = 800 s, 
and t = 10 3 s (from top to bottom) for xq = Mm. Temperature is drawn in units of 1 MK. The arrow below each panel represents 
the length of the velocity vector, expressed in units of 30 km s _1 . The color bar is common to all the panels. 
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Fig. 5. Same as Fig. 3, but for the case of oblique magnetic field (See Sect. 3.2). 
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Fig. 7. The SDO/AIA 304 A time sequence of a macrospicule propagation above the solar limb in the north-east polar coronal 
hole that clearly shows the upward motion of the plasma and its downward motion along the same path. 



